Phase separation in the Edwards model 
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The nature of charge transport within a correlated background medium can be described by 
spinless fermions coupled to bosons in the model introduced by Edwards. Combining numerical 
density matrix renormalization group and analytical projector-based renormalization methods we 
explore the ground-state phase diagram of the Edwards model in one dimension. Below a critical 
boson frequency any long-range order disappears and the system becomes metallic. If the charge 
carriers are coupled to slow quantum bosons the Tomonaga-Luttinger liquid is attractive and finally 
makes room for a phase separated state, just as in the t-J model. The phase boundary separating 
the repulsive from the attractive Tomonaga-Luttinger liquid is determined from long-wavelength 
charge correlations, whereas fermion segregation is indicated by a vanishing inverse compressibility. 
On approaching phase separation the photoemission spectra develop strong anomalies. 



I. PROBLEM 

The Edwards fcrmion-boson modcl^ constitutes a 
paradigmatic model for the theoretical description of 
quantum transport in solids. Charge transport normally 
takes place in the presence of some background medium^ 
Examples for such a "background" could be a spin-, 
charge- or orbital-ordered lattice^ but also a sequence 
of chemical side groups along the transport path, a de- 
formablc medium, or even a heat bath might be possible. 
In all these cases, the transfer of the charge carriers will 
be strongly influenced by fluctuations, which may exist 
in the background medium. The other way around, the 
properties of the background are quite often determined 
by the motion of the particle itself. 

Correlations inherent in such a complex interactive sys- 
tem are mimicked by a boson-affected hopping of spinless 
fcrmionic particles in the Edwards model. It reads: 

H = E Awl + M - A E( fe l + & i) + w ° E b t b i ■ 

(1) 

Every time a fermion fy' hops, it creates or absorbs a bo- 
son of energy ujq at the sites it leaves or enters. Such 
an excitation or de-excitation corresponds to a local "dis- 
tortion" of the background. Because of quantum fluctu- 
ations the distortions are able to relax (oc A). The physi- 
cally most interesting regimes in this setting are those of 
vanishing fermion density and of a half-filled band, rep- 
resenting doped Mott insulators, polaronic organics and 
charge-density- wave (CDW) systems^— with possible rel- 
evance to high-T c superconductors ) 10 ! 11 colossal mag- 
netoresistive manganites^— carbon nanotubeS ) 15 i 16 
graphen o 17 ' 18 and CDW transition metal complexes ji 9 .^ 
respectively. However, the Edwards model also reveals 
fascinating properties over the whole density range. 

On these grounds, the main goal of the present work is 
to pinpoint the ground-state phase diagram of the one- 
dimensional (ID) Edwards model as a function of the 



band filling n. Thereby we demonstrate that this model 
indeed captures a number of very interesting phenomena, 
including e.g. electronic phase separation (PS). To obtain 
reliable information about the ground-state and spectral 
properties of the model in the thermodynamic limit, we 
employ numerical pseudositc density-matrix renormaliza- 
tion group (DMRG) and dynamical DMRG (DDMRG) 
techniques (supplemented by a careful finite-size scaling 
procedure, for details see Refs.l9l. l2ll - |24| ). in combination 
with the analytical projective renormalization method 
(PRM),2^6 

So far the ID Edwards model has been solved ex- 
actly by numerical approaches for two casesi 22 i 24 ' 27 The 
first, in no way trivial case concerns just a single particle 
on an infinite lattice, where — depending on the model 
parameters — transport appears to be quasi-free, diffu- 
sive, or boson-assisted^^ 8 - When strong correlations de- 
velop at small A and large ujq, the background becomes 
"stiff" and the spinless particle's motion resembles those 
of a hole in an antifcrromagnctic spin background ) 29 ! 30 
as known from the t-J model. Interestingly, the Ed- 
wards model allows for so-called Trugman paths^I in an 
ID setting (Trugman paths usually unwind the string of 
misaligned spins a mobile hole leaves behind in two di- 
mensions). The second case is half- filling n = 1/2. Here 
a metal-insulator transition has been proven to exist: For 
small A and large loq the repulsive Tomonaga-Luttinger 
liquid (TLL) gives way to a CDW^ The CDW consti- 
tutes a few-boson state that typifies a Mott-Hubbard in- 
sulator rather than a Peierls state (normally established 
by the softening of a boson mode) . 



II. THEORETICAL RESULTS 

A. Ground-state properties 

The situation at finite density n (n ^= 1/2) is much less 
understood. By analogy with the t-J model one might 
expect that the system is metallic for < n < 1/2, 
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at least if the background is not too stiff. If so, the 
next question will be, whether the Edwards model might 
support the pairing of electrons in a certain parameter 
regime (A, u)q). Apparently a second electron, following 
the path of a first one, can take advantage of the back- 
ground excitations (bosons) left behind by the first elec- 
tron. This acts like an effective attractive interaction. 
If this attraction completely dominates the kinetic en- 
ergy the system might even phase-separate into particle- 
enriched and particle-depleted regions^ Since the Ed- 
wards model is not particle-hole symmetric, it is a moot 
point whether a corresponding hole pairing mechanism is 
at play also for 1/2 < n < 1. 

To address these problems, we analyze the charge cor- 
relations existing in the ground-state of the ID Edwards 
model. Firstly, we calculate the charge structure factor 
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for a system with TV fermions, L sites, and open boundary 
conditions (OBC). So the particle density is n = N/L 
and the momenta q = 2ns/ L with integers < s < L. 
The TLL charge exponent is proportional to the slope of 
S r ( q) in the long-wavelength limit q — > + (cf. Refs. [H 
andBl: 



K„ = 7r lim 



S c (q) 



(3) 



Then, K p > 1 (K p < 1) characterizes an attrac- 
tive (repulsive) TLL for our spinless fermion transport 
model ((T|), and K p = 1/2 will define a metal- insulator 
transition point at n = 1/2M- That the finite-size scaling 
of S c (q) and K p works well has been demonstrated for 
the half-filled band case^ Secondly, in order to comprise 
potential PS, we determine the finite-size equivalent of 
the charge compressibility^ 
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(4) 



with Eq(N) being the ground-state energy for N elec- 
trons on L sites. An infinite compressibility signals the 
formation of a PS state. 

Fig. [U demonstrates that all this works in practice. Ex- 
emplarily choosing n = 1/8 and A = 0.2, we show how 
the TLL exponent K p and the inverse compressibility 
k" 1 scale with increasing system size at various ujq. The 
transition point between a repulsive TLL — obviously re- 
alized at luq > wo, c i when excitations of the background 
are energetically rather costly — and an attractive TLL at 
smaller ojq, can be read off from the inset [depicting the 
extrapolations K p (L -» oo)] tobewo,ci(A = 0.2) ~ 1.118. 
Reducing luq further, in the attractive TLL phase, a dra- 
matic increase in K p is observed at wo,c2 — 0.6. Our in- 
verse compressibility data indicate that at this point the 
attraction among the particles becomes so strong that 
the system shows PS, i.e., n^ 1 = for u>q < ujq C 2- 
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Figure 1. Finite-size scaling of the TLL parameter K p (left 
panel) and inverse compressibility k~ (right panel) for the 
ID Edwards model with A = 0.2 at n = 1/8. Insets show 
the variation with uio of the L — > oo extrapolated values. 
Results are obtained by DMRG for a lattice with L sites and 
OBC. In all DMRG runs we take into account up to n;, = 
4 bosonic pseudosites and determine rib by the requirement 
that the local boson density of the last pseudosite is less than 
1CP 6 for all j. Using selectively n;, = 5, we controlled that 
convergence is truly achieved. Furthermore we keep up to 
m = 1200 density-matrix eigenstates in the renormalization 
steps to ensure that the discarded weight is smaller than 10~ 8 . 



Proceeding in the same manner for different values of n 
and A respectively uio, we can map out the phase diagram 
of the ID Edwards model. The outcome is displayed in 
Fig. [2] Let us first consider the case of a not too small 
boson relaxation (oc A), which ensures that the system is 
metallic for large ojq in the whole density regime. Then, 
as the upper panel of Fig. [5] shows, depending on n we 
find up to three different regimes: For small and large 
particle densities a repulsive TLL, an attractive TLL and 
a PS state appear in sequence as the energy of the bosons 
is lowered. In contrast, around half- filling only the repul- 
sive TLL exists. At this point we would like to remind 
the readers that at half-filling, for smaller values of A, 
(A < A c ), a CDW is formed for wo > uo, c ^ Our unbi- 
ased DMRG calculations give no evidence for any other 
phases. Note that the behavior in the low-density regime 
is consistent with what is found for the ID t-J mode l 36 ' 37 , 
where the holes correspond to the spinless fermions in the 
Edwards model. As mentioned before, the phase diagram 
is not symmetric with respect to n = 1/2. This is because 
the hopping of a hole (missing electron) to a neighboring 
site creates in the Edwards model a boson at the arrival 
site and not on the departure site as in the motion of 
a single particle. Since this boson can be destroyed im- 
mediately when the hole makes a further hop, there is 
no string effect and a few holes propagate more easily 
than a few electrons. Hence the background at fixed ujq 
and A appears to be less stiff for n < 1 than for n > 
and, as a result, the boundary between the repulsive and 
attractive TLL is shifted to larger values of CJo >c i- Con- 
trariwise, since carriers will be less mobile in a state with 
charge separation, the boundary between the attractive 
TLL and the PS state wo. C 2 is shifted to smaller values if 
one compares the corresponding results at high and low 
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carrier density. 

We next consider a fixed particle density n — 1/8 and 
track the phase boundaries in the \-uiq plane; see the 
lower panel of Fig. [5J The repulsive TLL established for 
large ljq is strongly renormalized if the ability of the back- 
ground to relax is low. For example, we find K p = 0.527 
for lj = 2 at A = 0.01. Lowering uj for such a stiff 
background, the transition to the PS state happens al- 
most instantaneously with a narrow intervening attrac- 
tive TLL state; at A = 0.01 we have Wrj,ci — 0.747 and 
o;o ;C 2 — 0.672. If we fix, on the other hand, ljq = 2 and 
increase A, we observe a strong enhancement of K p in 
the interval < A < 0.3 (K p = 0.767, 0.879, 0.932 for 
A = 0.1, 0.2, 0.3, respectively), followed by a very gradual 
increase until the transition to the attractive TLL takes 
place at about A c i — 1. Obviously, the region where the 
attractive TLL (PS state) constitutes the ground state 
expands (shrinks) as A gets larger, which can be traced 
back to the subtle competition between kinetic energy 
and charge segregation effects. 

To gain deeper insight into the different mechanisms 
at play, we analyze the variation of the local fcrmion 
[rij = {fjfj)] and boson [n* = (b^bj)] densities in real 
space. In Fig.[3]we show both density profiles for n = 1/8 
and characteristic values of A and wq, implementing re- 
pulsive TLL, attractive TLL, and PS states. In the 
metallic regime, the OBC lead to Friedel oscillations in 
the particle density. These oscillations are especially pro- 
nounced for the repulsive TLL [panel (a)], and they will 
be even stronger for ujq = 2 (not shown) , where the num- 
ber and fluctuations of the bosons is reduced. Note that 
the Friedel oscillations caused by the OBC will be alge- 
braically reduced, if we move for very large systems to- 
wards the central part of the chain. If we enter the attrac- 
tive TLL regime by increasing A at fixed u>o the Friedel 
oscillation will be smeared out [panel (b)]. Thereby, the 
number of oscillations stays constant, which means that a 
pairing of electrons does not occur. This is an important 
difference from the spinful ID t-J model, where a recent 
DMRG study reveals that the number of Friedel oscilla- 
tions is halved (by increasing J), corresponding to half 
the number of particles, which indicates pairing 4£ Next 
lowering luq to a point where PS sets in, the particle den- 
sity oscillations vanish, see panel (c). There, just above 
the PS boundary, no evidence is found for the clustering 
of multiple particles in several groups. This agrees with 
the findings for the ID t-J models Once we are deep 
inside the PS region, a density distribution results, where 
a single island of particles at the central part of the sys- 
tem appears, leaving a sizable number of almost empty 
sites at both ends of the chain [see panel (d)]. Looking 
at the bosonic degrees of freedom, we see that the strong 
attraction among these particles is mediated by a boson 
cloud covering the electron-rich region. As a result the 
kinetic fluctuations will be strongly quenched. 

It is all but impossible to comprise by DMRG the large 
number of bosons and the strong bosonic fluctuations, 
which appear at still smaller loq in the Edwards model, 
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Figure 2. Ground-state phase diagram of the ID Edwards 
model, showing repulsive TLL (K p < 1), attractive TLL 
(K p > 1), and PS (k _1 = 0) regions. The phase boundaries 
were obtained by DMRG (filled symbols) in the course of a 
careful finite-size scaling analysis and for the infinite system 
directly by PRM (open symbols). The upper panel displays 
the phase diagram as a function of n, varying loq at fixed 
A = 0.2; the lower panel gives the phase diagram in the X-cjq 
plane for fixed density n — 1/8. 



simply because the dimension of the Hilbcrt space in- 
creases dramatically. To access the low-wo regime and 
reconfirm the DMRG phase boundaries, obtained for the 
ID Edwards model at larger values of wo, we employ 
the analytical PRM approach^ The basic idea of the 
PRM is to construct — performing a sequence of discrete 
unitary transformations and eliminating all transitions 
with energy larger than a given cutoff energy — an effec- 
tive non-interacting Hamiltonian H with renormalized 
parameters (in the limit of vanishing cutoff energy). For 
the metallic state of the Edwards model, in momentum 
space, it takes the form 



E, 



(5) 



where the renormalization equations for ik and w g have 
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Figure 3. Local particle (filled circles) and boson (open 
squares) densities for the ID Edwards model with n = 1/8, 
calculated for a 48-site system with OBC at various A and luq. 



been derived in Ref. [26| (for the half- filled band case). 
In order to fix the boundary between metallic and PS 
states, we (i) solve the renormalization equations at a 
given Ep in the TLL phase, (ii) determine the corre- 
sponding fermion density n, and (iii) slightly vary Ep 
(to get closer to the PS instability), determine n, and re- 
peat the whole procedure self-consistently until the func- 
tional dependence of n on Ep is established. A vanish- 
ing inverse compressibility can then be read off from a 
sudden jump of n under a tiny variation of Ep or, the 
other way around, if a plateau in the Ep versus n plot 
appears. The transition points, determined in the man- 
ner described, were inserted in Fig. [2] to complete the 
ground-state phase diagram in the n-u)Q plane. Whenever 
points can be compared we find a remarkable agreement 
between our DMRG and PRM data. The deviation of 
the PRM data at larger particle densities n > 0.75 re- 
sults from uncertainties in fixing the jump of n under 
Ep variation. Here we should trust in the DMRG phase 
boundary. 



B. Spectral properties 

Let us finally discuss the single-particle spectrum of 
the ID Edwards model. The single particle excitations 
associated with the injection (+) or emission (— ) of an 
electron with wave vector k, 



A ± (k,tv) 



E 



\(i>t\fk\M a ^Tu> ± ), (6) 



can be computed by DDMRG^^ where f£ = fl, = 
f k , \ipo) is the ground state of a L-site system in the TV- 
particle sector, and \4>m) denote the m-th excited states 
in the (TV ± l)-particle sectors with excitation energies 
w± = E± -E . 

Within PRM we find for the photoemission part 

A~(k,ui) = 5%, hf, 6(cj — £fe) 

+ ^lk, q n b q n f k _ q S(u}-Q q -e k - q ), (7) 

where n£ (h q ) are the fermion (boson) occupation num- 
bers in momentum space calculated with the renormal- 
ized Hamiltonian H. The coefficients a|, P\ q , and 
7^ follow from renormalization equations^ Taking the 
corresponding expression for A + (k,uj), the sum rule 
JT~ duj[A+(k, uj) + A~(k,u)] = 1 is fulfilled. 

Figure 3] gives the combined photoemission spectrum, 
A(k,w) = A + (k,w)+A~(k,u}), as obtained by DDMRG 
for quasimomenta k. In the repulsive TLL regime (left 
panel, loq = 2) we find a rather coherent signal with 
comparable spectral weight for all k values k < 7r/2. For 
larger k, excitations with at least one additional tJo-boson 
involved become important (recall that both initial TV- 
particle- and target (TV ± l)-particle states are multibo- 
son states with a momentum being the total momentum 
of electrons and bosons). The spectrum in the attractive 
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Figure 4. Line-shape of the single-particle spectral function 
A(k,u>) for the ID Edwards model with A = 0.2 at n = 1/8. 
Results are obtained by DDMRG for an 16-site chain with 
quasimomenta k = %s/(L + 1), using a broadening r\ — 0.2. 
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Figure 5. PRM [inverse] photoemisson spectrum [A + (fc,cj)] 
A~(k,cj) for the Edwards model with A = 0.2, n = 3/10. The 
bold line marks the signal at the Fermi momentum JtF = O.Sn. 

TLL phase (right panel, ojq = 0.7) shows a sharp absorp- 
tion signal in the vicinity of kp only. Here the (inverse) 
photocmission spectrum for k < kp (k > kp) exhibits a 
few absorption maxima at multiples of the boson energy. 
Obviously, due to the smaller u , here the dynamics of 
the system becomes dominated by bosonic fluctuations. 

The (inverse) photoemission close to the transition to 
the PS state is depicted Fig. as calculated by PRM 
for wo = 0.5. Here the almost dispersionless signal for 
k < 0.2127T (k > 0.7887T) is a precursor of the PS instabil- 
ity. As we noted above, Ep docs not show any variation 
with n if n is smaller (larger) than the lower (upper) crit- 



ical density n C 2 for PS. Since in the PRM derivation of 
A^(k,ui) terms with two (and more) bosonic creation or 
annihilation operators were neglected,— the photoemis- 
sion spectrum of Fig. [5] does not feature the multiboson 
related signatures found within a DDMRG treatment. 

III. CONCLUSIONS 

To summarize, the combination of analytical (PRM) 
and numerical (DMRG) approaches permits the precise 
determination of the ground-state phase diagram of the 
ID Edwards model in the whole parameter regime. In 
the low- and high-density regions, the attraction between 
the particles mediated by the bosonic degrees of freedom 
representing the background medium might become so 
strong that electronic PS sets in. In the remaining re- 
gion the system realizes a TLL or, at n = 1/2, possibly 
even a truly long-range ordered CDW state. Depending 
on the properties of the background medium, the TLL 
might be attractive or repulsive. The richness of the 
phase diagram is remarkable. The model captures im- 
portant features of Holstein, t-J, Hubbard, and Falicov- 
Kimball type models^— >2£ Since the Edwards model is 
one of the simplest models for studying transport in low- 
dimensional systems, an inspection of its predictions with 
ultracold fermion/boson quantum gases^S as, e.g., car- 
ried out for the ID t-J model^i would be of great inter- 
est. 
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